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Abstract: Unusual emission of light, called the unconventional Smith- 
Purcell radiation (uSPR) in this paper, was demonstrated from an electron 
traveling near a finite photonic crystal (PhC) at an ultra-relativistic velocity. 
This phenomenon is not related to the accepted mechanism of the conven- 
tional SPR and arises because the evanescent light from the electron has 
such a small decay constant in the ultra-relativistic regime that it works 
practically as a plane-wave probe entering the PhC from one end. We 
analyze the dependence of the SPR spectrum on the velocity of electron 
and on the parity of excited photonic bands and show, for PhCs made up of 
a finite number of cylinders, that uSPR probes the photonic band structure 
very faithfully. 
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1. Introduction 

A traveling charged particle induces coherent radiation when it passes near a periodic dielectric 
structure along the direction of its spatial periodicity. This radiation, called Smith-Purcell radi- 
ation (SPR) [1,2], can be a novel radiation source with several remarkable properties. The most 
important one is the scalability of the output frequency; the threshold frequency below which 
SPR is kinetically impossible varies in inverse proportion to the magnitude of the period. In ad- 
dition, the SPR is characterized by the presence of resonances at a series of frequencies, which 
again vary in inverse proportion to the period. Owing to these properties, it has been recognized 
that SPR can be a basic mechanism for a compact free-electron laser [3]. 

Since its first observation, SPR has been studied using mainly metallic diffraction gratings 
of one-dimensional periodicity [4, 5, 6, 7, 8]. In most theoretical analyses made so far, the 
gratings have been treated as perfect conductors to simplify the treatment of the periodic light 
scattering [9, 10, 11, 12]. In these systems, Wood's anomaly [13] in the optical density of states 
(ODOS) is responsible for the enhanced signals of SPR, and thus the relevant frequencies of 
the resonances are determined in a straightforward manner using simple kinetics. 



Recently, both theoretical [14, 15, 16, 17, 18, 19, 20] and experimental [21] SPR results 
have been reported for photonic crystals (PhCs) used in place of metal gratings. It was found 
that PhCs induced highly coherent SPR because of their multidimensional periodicity in their 
dielectric functions. The SPR spectrum consists of point-like signals as a function of frequency, 
which show up each time the evanescent light from the electron excites a photonic band (PhB) 
mode of high quality factor. SPR from a PhC is versatile, because PhCs generally have various 
parameters, which can now be reliably designed and changed. 

However, when an electron beam of ultra-relativistic velocity was used in combination with 
a PhC, which was finite in the direction of the electron trajectory, unexpected phenomena that 
contradicted the conventional understanding of the SPR were experimentally observed [22]. 
Such phenomena have not been observed in the gratings of nearly perfect conductors and are 
expected to be absent even in PhCs when an electron beam of slower velocity is used. This 
SPR, called unconventional SPR (uSPR) in this paper, is quite distinct in many ways from the 
conventional SPR (cSPR) and thus is easily identified; most importantly, the uSPR spectrum 
sweeps the entire region of frequency-momentum phase space, in contrast to the cSPR which 
carries information of the phase space only along the shifted v lines, to be defined later. In the 
phase space, the uSPR is characterized by peculiar resonances arising along curves, which are 
related more or less to the dispersion relations of PhBs, not to the shifted v lines, with relatively 
little intensity variation. Therefore, SPR in the ultra-relativistic regime of the beam velocity is 
potentially useful both as a monochromatic light source and as a probe to investigate the PhB 
structure. We should note here that the Cherenkov effect also can be used to probe the PhB 
structure with the angle-resolved electron energy loss spectroscopy [23, 24, 25]. 

The conjecture inferred by the experiment, which this paper seeks to verify, was that the SPR 
consists of the cSPR and uSPR in the relativistic regime of electron velocity and that the uSPR 
is expected to disappear gradually with decreasing velocity, to leave solely the cSPR component 
for velocities typically less than a few hundred keV. 

Very recently, Kesar et al. [26] reported a systematic discrepancy between the calculated 
SPR spectrum in a finite-size grating and that of the conventional theory assuming infinite size. 
Since they focused on the diffraction grating of a perfect conductor, their discrepancy is not 
directly related to ours, which arises in systems involving PhCs with a finite dielectric constant. 
As for this point, it is noteworthy that a rigorous theory was developed for finite-size grating of 
infinitely-thin metallic plates [2]. 

This paper presents a comprehensive theoretical analysis of the uSPR for the PhCs composed 
of a finite number of cylinders. We use the multiple-scattering method of theoretical treatment 
[19], which explicitly takes into account the finiteness of the total number of cylinders and 
treats exactly the multiple Mie scattering among them. We shall investigate the properties of 
the uSPR and the interplay between u and c SPRs by changing various parameters, such as 
electron velocity, dielectric constant of PhCs, length of PhCs, thickness of PhCs and angle of 
SPR emission. Most of this paper will focus on the PhCs of dielectric cylinders of circular cross 
section. The remarkable difference of the present work from the previous theoretical works 
involving PhCs lies in the fact that we are dealing with the finiteness of the length of PhC in the 
direction of the electron beam. 

This paper is organized as follows. In Sec. II, we present the kinetics for both cSPR and uSPR 
and discuss what plays a key role in inducing uSPR. Section III is devoted to the comparison 
of the spectrum of uSPR with that of cSPR. In Sec. IV, we present detailed analyses of uSPR 
by changing various parameters. Finally, we summarize the results in Sec. V. 
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Fig. 1. Schematic illustration of the system under study. An electron travels with constant 
velocity v and impact parameter b below the one-dimensional periodic array of dielectric or 
metallic cylinders. The cylinders are arrayed periodically with axes in the z direction with 
radius r, dielectric constant e and the lattice constant a. The trajectory of the electron is 
parallel to the direction of the periodicity. Angles 6 and <p are the polar angles of the SPR 
signals. This PhC has a mirror plane indicated by dotted line in the right panel. 



2. Kinetics of conventional and unconventional SPRs 

In the following discussion, we focus on a PhC composed of infinitely long cylinders. The 
system under study is schematically illustrated in Fig. 1. Cylinders are arrayed periodically 
(lattice constant: a) in the x direction with cylinder axes in the z direction. An electron travels 
near the PhC in a trajectory parallel to the x axis with velocity v and impact parameter b. 
We obtain the SPR spectrum from this system as a sum of the plane-wave signals generated 
by the scattering of the evanescent light emitted by the electron. The whole process of multiple 
scattering among a finite number of cylinders is dealt with compactly by the multiple-scattering 
theory of radiation using the vector cylindrical waves as a basis of representation [19]. 

Let us briefly summarize the kinetics in the theory of cSPR. In the theory of cSPR, a finite 
periodic structure is simulated by a periodic structure of infinite length in the x direction. A 
traveling electron accompanies the radiation field that is a superposition of evanescent waves 
with respect to frequency CO and wave number k z in the z direction [16]. The wave vector of 
each evanescent wave is given by 

where F is purely imaginary because v <c. The imaginary part |T| determines the spatial decay 
of the evanescent wave incident on the PhC. In what follows, it is important to remember 
the feature of k z = that, in the ultimate limit v — > c, |T| tends to zero. Since CO and k z are 
conserved quantities in the geometry of Fig. 1, the evanescent waves with different CO and k z 
are independent in the whole scattering process. Thus, the incident light of CO and k z leaves the 
PhC, after being scattered, with the same CO and k z . Therefore, the SPR signals observed in the 
xy plane may be analyzed by setting k z = everywhere. Since we are now dealing with a perfect 
periodicity extending from — °° to °°, we obtain the SPR signal in the form of Bragg-scattered 
waves summed over the diffraction channels. The channels are specified by the wave vector 
defined by 

K± = (^-h,±r h ,k z ), r^^) 2 -^-*) 2 -**. (2) 

Here, h = Inn /a (« : integer) is a reciprocal lattice point of the PhC in the x direction. Before 
the scattering by the PhC, CO and k x , the x component of the wavevector of light, satisfy the 
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Fig. 2. Schematic illustration of the input evanescent light yielding conventional (left 
panel) and unconventional (right panel) SPRs. The conventional SPR is produced when 
the evanescent wave has an appreciable decay constant. The incident light enters the PhC 
from below. The unconventional case arises when the incident evanescent light has a neg- 
ligible decay constant and is regarded as a plane wave entering the PhC from its left edge. 
In this case, the evanescent wave is almost symmetric with respect to the mirror plane, 
inducing the even-parity selection rule in the PhB excitation. 

relation k x = co /v. The line CO — vk x , called the v line in this paper, lies outside the light cone in 
the phase space (k x , co). After the scattering, the x component of the light of channel h becomes 
k x = co/v — h. The shifted v line defined by this equation is inside the light cone in a certain 
frequency range. In that frequency range we can detect the SPR signal in this channel at a 
far-field observation point. The propagating direction of the SPR signal of co is given by 



in the polar coordinates defined in Fig. 1. 

The inside region of the light cone is the leaky region of PhBs, and, accordingly, the ODOS 
is nontrivial there. Actually, ODOS has a sequence of peaks of finite width in the (k x ,k z ,Co) 
space. The peak position determines the dispersion relations co = co„(k x ,k z ) of the quasi-guided 
PhB modes. Imagine temporarily k. = 0, for brevity. The presence of the modes in the (k x ,a>) 
space significantly affects the SPR spectrum by causing a sharp resonance when the dispersion 
curves of PhBs intersect the shifted v lines. The resonance becomes sharper as the quality factor 
of the relevant PhB modes increases [16, 17]. Therefore, the SPR from a PhC can have very 
high quality, which is intriguing as a new possibility of PhC. 

In the above argument, we assumed the conservation of k x , with an Umklapp allowance 
taken into account. In an actual PhC with a finite number N of cylinders, the periodicity or the 
translational invariance of the whole system is lost at the sample edges. One way to take account 
of the finiteness of N is to treat k x as defined only approximately with a width of the order of 
Ak x ~ 2n/(Na) [2]. In this approach, the shifted v lines are considered to have the finite width, 
and the PhB dispersion relation will become detectable within this allowance centering on the 
shifted v lines of open channels. In reality, however, the uSPR signals appear in the phase space 
(k x , co) with a much larger distribution than this straightforward 1 /N blurring [22]. 

Let us now consider an ultra-relativistic velocity v ~ c. We continue to confine ourselves 
to the measurement within the xy plane, i.e., k z = 0. In this case, T is almost zero, and the 
evanescent wave incident on the PhC may be regarded practically as a plane wave with its 
wavevector directed in the x-direction. Therefore, the light-scattering problem is quite similar 
to that of the light transmission and reflection in the x direction through the periodic array of 
cylinders, as depicted in Fig. 2. In this situation, it is obvious that the sample edges play a 
crucial role. In particular, we know that k x is no longer a good quantum number, no matter how 
large the total number of the cylinders may be. As in an ordinary light-transmission experiment 
involving PhCs, we expect the incident wave of frequency CO to excite the PhB modes at the 
crossing points between the line of constant co and the dispersion curves CO = (0„(k x ,Q) of PhBs 
in the (k x , co) plane. The point is that the wave vector of the PhB thus determined is not related 




(cos , sin cos (j) , sin 6 sin <p ) , 



(3) 



to the value of k x on the v line of the incident light. Therefore, we should expect SPR signals 
over the entire (k x ,(o) space, not necessarily restricted along the shifted v lines. The signals 
expected off the shifted v lines characterize the uSPR. Also, we can expect that the transmitted 
SPR obtained in the side opposite to the trajectory is almost identical to the reflected SPR of 
the trajectory side. Finally, similar to the ordinary setup of a plane wave transmitting through 
a two- and three-dimensional PhC [27], a selection rule must exist for the symmetry of the 
PhB modes to be excited. The PhC in our problem has a symmetry with respect to the mirror 
reflection y — > — y (y = is the plane bisecting the PhC), and in the ultra-relativistic regime the 
incident light is of y independent. Hence, solely the PhB modes of even mirror-symmetry are 
expected to participate in the resonant light scattering. The even-parity selection rule will thus 
characterize uSPR spectra. 

According to this scenario, a PhB mode manifests itself in the uSPR. Therefore, the uSPR 
will have a rather broad band as a function of frequency. This is in contrast to the cSPR, in 
which excited PhB modes give rise to sharp resonance peaks only along the shifted v lines. In 
addition, since v ~ c, the shifted v lines coincide with the threshold lines for the opening of a 
new Bragg diffraction channel. The channel opening often leaves a singular trace due to Wood's 
anomaly in the line shape of wave scattering. Thus, on our shifted v lines, Wood's anomaly will 
occur, together with the resonance peaks of the cSPR associated with the PhB excitation. In 
this way, the spectra of the c and u SPRs reveal a quite rich structure when the electron is 
ultra-relativistic. 

So far, we have concentrate ourselves on the electron traveling parallel to the x direction, that 
is, the direction of the periodicity of the PhC under consideration. If the velocity vector v of the 
electron is given by 

v = (v v ,0,v,) = v(cos a, 0, sin oc), (4) 

the kinetics of SPR changes accordingly. In particular, the dominant component of wave num- 
ber k z depends on frequency and is given by k z = v-o/v 2 . Within the conventional theory the 
SPR acquires a significant enhancement when the following three conditions are fulfilled: 

a> = v x {k x - h) + v z k z , k z = -^(o, a> = (0„(k X7 k z ), (5) 

v 

where the first equation defines the shifted v line at nonzero v z . As in the case of v z = (a = 0°), 
this scenario of the SPR is insufficient for an ultra-relativistic electron. In this case the evanes- 
cent wave accompanied by the electron can be effectively treated as a plane wave propagating 
parallel to v. This highlights the role of the sample edge of the finite-size PhC, namely, the bro- 
ken translational invariance, and the photonic band modes on the entire plane of k z = (v z /v 2 )(0 
are excited, not necessary restricted on the shifted v lines. In the limiting case of vanishing 
v x (a = 90°), the electron travels parallel to the cylindrical axis. There, no propagating radi- 
ation is generated from the PhC, as far as the cylinders have infinite length in the z direction. 
This is due to the perfect translational invariance along the axis. In actual PhC, however, this 
translational invariance is broken, yielding a sort of the diffraction radiation. Thus, as a varies 
from 0° to 90°, the conventional theory of the SPR, which assumes the translational invariance 
both in the x and z direction, predicts a gradual disappearance of the SPR. On the other hand, 
the uSPR gives a novel radiation irrespective of a, in which the broken translational invariance 
in the x direction is highlighted at small a, and that in the z direction is highlighted around 
a = 90°. 

To summarize, the two SPR spectra, c and u SPRs, coexist in the ultra-relativistic regime, 
when a PhC sample is used. To analyze the experimental signals based on the knowledge on 
the band structure of photons, the length of the PhC or the total number of cylinders must be 
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Fig. 3. PhB structure of the monolayer of low-index cylinders in contact (e = 2.05). The 
modes of TE polarization of k z = are plotted as a function of k x . The PhB modes are 
classified according to the parity with respect to the mirror plane bisecting the monolayer. 
The light line CO = ±ck x is indicated by thick solid lines, the shifted v lines of v = 0.5c by 
thin solid lines, and those of v = 0.99999c by dashed lines. The horizontal arrows (six for 
v = 0.99999c and four for v = 0.5c) are drawn at the intersections between the shifted v 
lines and the PhB dispersion curves. They correspond to those of Figs. 4 and 7. 

finite but large enough. For a finite system to have a band structure comparable to that obtained 
for an infinite system, the periods of, N > 8 will be enough according to our experience. In 
contrast to the value of N in the x direction, however, we are considering a system having 
small size in the y direction, such as a PhC made of a monolayer or stacked layers of several 
monolayers. The finite size in the y direction needs to be considered explicitly to obtain the 
band structure of our PhCs. 

3. Typical example of conventional and unconventional SPRs 

Before making a detailed study of the uSPR properties, we briefly compare the cSPR spectrum 
with the uSPR spectrum, using the numerical results for a test system. We adopt a PhC of a 
monolayer of periodic array of low-index cylinders (dielectric constant e = 2.05). For a radius- 
to-periodicity ratio r/a = 0.5, Fig. 3 depicts the band structure of the monolayer of an infinite 
number of cylinders. The band structure inside the light cone was obtained by plotting the 
peak frequencies of the ODOS, which were calculated as a function of k x and k z [28]. The 
band structure outside the light cone (that of the true-guided modes with infinite lifetime) was 
obtained from the position of the poles of the S matrix, which are found on the real axis of the 
complex co plane. In Fig. 3, k z — is assumed, so that the PhB modes are decomposed into 
purely transverse-electric (TE) and transverse-magnetic (TM) modes. Only the band structure 
of the TE modes is presented, because the incident evanescent wave is TE -polarized at k z = 0. 
The PhB modes are further classified by parity with respect to the mirror plane y = 0. In Fig. 
3 the even (odd) parity modes are indicated by red (green) circles. We should note that PhB 
modes having their dispersion curves disconnected in Fig. 3 are the ones obtained from the 
ODOS peaks, which are often too broad to identify the peak position. 

Let us first consider the cSPR spectra obtained by the conventional theory based on the 
assumption of perfect periodicity in the x direction. We used the parameters v = 0.99999c, 
<p = 180° and b = 3.33a and assumed that the radiation was observed in the xy plane (k z = 0), 
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Fig. 4. Reflected intensity of cSPR on the shifted v lines of v = 0.99999c for £ = 2.05 and 
b = 3.33a. Perfect periodicity ranging from x = — °° to x = °° of the monolayer cylinders is 
assumed. The arrows are drawn at the peak positions and agree with those of Fig. 3, which 
were assigned to the crossing points between the shifted v line and the band structure. 

as actually encountered in the millimeter-wave SPR experiments carried out recently [21, 22]. 
Since the periodicity is perfect, the SPR spectrum appears strictly on the shifted v lines, which 
are almost parallel to the light line CO = ck x . Figure 4 presents the reflected cSPR spectra along 
the shifted v line of h = 1 and 2 (in units of 2n/a). The peaks of the cSPR spectrum arise at 
the frequency where the shifted v lines k x = co/v — h intersect the PhB structure given in Fig. 
3. Several arrows are drawn at the peak positions in Fig. 4 and, to identify each of the peaks, 
horizontal arrows are added in Fig. 3 at the corresponding positions in phase space. Comparing 
these two figures, we see that the peak lowest in frequency arises from the excitation of an odd- 
parity PhB mode. Thus, the even selection-rule for the parity of the excited PhB modes does 
not hold for cSPR, though it somewhat affects their spectral shapes. 

Now we turn to the uSPR spectrum. The SPR spectrum, with the finiteness of N considered 
explicitly, is obtained over the entire (k x , co) space by summing all the amplitudes of the multi- 
ply scattered light from the N cylinders [19]. The numerical result for N = 21 is given in Fig. 5, 
for the same parameters of r/a and e as used in Fig. 4. The angle 9- and frequency 0)-resolved 
reflected SPR intensity is mapped onto the (k x , co) plane through the relation k x = {co/c) cos 9. 
To be precise, |/ M (0)| 2 with -% < 9 < defined in Eq. (33) of Ref. [19] was plotted by using 
the above relation. 

We observe that the peaks of the SPR intensity are found along 

(A) the shifted v lines, 

(B) the curves whose slopes are positive and less than 1, 

(C) the curves whose slopes are negative, 

(D) the forward light-line (to = ck x ), and 

(E) the fiat lines terminated on the backward light-line (co = —ck x ). 

For comparison, we superimpose the even-parity PhB structure of TE modes on the intensity 
map and present the result in Fig. 6. From Fig. 5, we find that type (A) peaks are rather broad 
along the shifted v lines as a function of frequency, as compared to the cSPR spectrum shown 
in Fig. 4. A good agreement between the PhB dispersion curves and the high intensity positions 
of the uSPR demonstrates that the peaks of type (B) are attributed to the quasi-guided PhBs of 
even parity. Also, no evidence of the excitation of the odd-parity modes is found in Figs. 5 and 
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Fig. 5. Reflected SPR intensity map from a finite (N = 21) monolayer of the contact cylin- 
ders with a low index dielectric constant. The signals appearing off the shifted v lines 
characterize the uSPR. Except for A', the same parameters as in Fig. 4 were used. 




Fig. 6. Figure 5 of the intensity overlaid with the PhB structure of the even parity (indicated 
by red circles). The shifted v lines are indicated by solid lines. 
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Fig. 7. The reflected cSPR intensity from the infinite monolayer (N = °°) of the low-index 
contact cylinders. The parameters v = 0.5c and b = 0.2a were assumed for the electron 
beam. Four arrows are assigned to the peak positions and coincide with those of Fig. 3. 

6. Thus, the even-parity selection rule of uSPR, predicted in the previous section, indeed holds. 
The peaks of type (C) are found, for instance, around (k x a/2n,COa/2nc) = (—0.75,1.25), 
whose origin is also attributed to excitation of the even-parity PhB modes. The intensity of 
these peaks is, however, small, as compared to that caused by the PhB modes with positive 
slope. This seems reasonable, considering that the PhBs of the negative group velocity will 
have suppressed excitation probability at the left edge of the PhC. The peaks of type (D) are 
inevitable in finite systems; the incident wave induces the quasi-guided waves in the monolayer, 
which exit from the right edge, causing a forward-oriented diffraction there. As a consequence, 
broad peaks of the SPR near the forward light-line emerge. Note that the intensity of the type 
(D) signal oscillates as a function of frequency. This is a Fabry-Perot oscillation of the signal 
intensity with period decreasing with increasing N. Finally, the signals of type (E) are caused 
by the presence of pseudo gaps in the PhB structure. To see this, we have only to note in Fig. 
6 that the flat streaks of type (E) signals appear at the gap positions. When the frequency of 
the incident wave lies in a pseudo gap of the monolayer, the incident wave cannot penetrate 
deep inside the PhC and is scattered out of the PhC as a type (E) signal, with a certain angular 
distribution centering on the backward direction. Note that the intense streaks are not related to 
the band gaps of odd parity. 

The above features of the uSPR in finite monolayers will remain unchanged even for a semi- 
infinite monolayer, which is made of cylinders ofN = °° but bounded at one end, because what 
matters in the above discussion is the presence of the left edge of PhC as an entrance surface of 
a wave propagating in the x direction. 

Next, we consider SPR spectrum obtained from a slower electron in a non-relativistic regime. 
The calculation is made for v = 0.5c, which is a typical value for the electron velocity used in 
scanning electron microscopes. The parameters except v and b are the same. As above, we 
compare two spectra of N = °° and N — 2\. 

First, the reflected SPR spectrum for N = °° is given in Fig. 7. The spectra reveal a marked 
resonance at a>a/2nc = 0.621. The line shape of the resonance is asymmetric as a function of 
frequency. As indicated by arrows, each agreeing precisely with those given to the shifted v line 
of Fig. 3, the cSPR peaks all appear exactly at the intersection points of the shifted v line of 
v = 0.5c with the PhB dispersion curves. 
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Fig. 8. Reflected SPR spectrum from a finite (jV = 21) monolayer of low-index cylinders 
in contact. The intensity profile is overlaid with the corresponding PhB structure of N = °°. 
The PhB modes with even (odd) parity are indicated by red (green) circles. The parameters 
v = 0.5c and b = 0.2a were assumed for the electron beam. 

The reflected SPR spectrum for N = 21 is given in Fig. 8, with the superposition of the PhB 
structure (of N = °°). We see at once that high intensity SPR appears only on the shifted v lines, 
although very weak structures reminiscent of the finiteness of our PhC are still seen off the 
shifted v line. This is in clear contrast to the ultra-relativistic spectra, where marked signals of 
uSPR existed definitely off the shifted v lines. The signals on the shifted v line of h = 1 have a 
resonance peak at (k x a/2n, coa/2nc) = (0.24,0.62). This frequency is almost identical to that 
of the resonance obtained for N = °° shown in Fig. 7. Also, we can perceive the asymmetry 
of the line shape along the shifted v line, as in the cSPR spectrum of N = °°. Therefore, we 
may conclude that, for slower velocities such as v = 0.5c, the SPR of the finite PhC can be 
understood sufficiently well using the theory of cSPR, based on the assumption N = °°. The 
uSPR signals are suppressed as follows. The light of v = 0.5c is literally evanescent with an 
appreciable decay constant |r|, so that, while passing through the PhC in the +y direction, 
the incident light decays much and sees only the surface region of cylinders. Accordingly, the 
picture of a plane wave with wavevector in the x direction no longer holds and the conventional 
theory of SPR covers all the features. 

4. Properties of the unconventional SPR 

This section presents the properties of uSPR in detail by changing various parameters. As ex- 
plained in the previous sections, the broken translational invariance due to finite number of 
cylinders (AO is crucial in the uSPR. Taking account that the uSPR must vanish in the system 
of the perfect translational invariance, it is interesting to investigate the A^-dependence of the 
uSPR in detail. The number of stacking layers (N/) is also an important factor because ODOS 
and thus the PhB structure depends crucially on N/. Dielectric constant e and radius r of the 
cylinders are other factors that significantly influence the PhB structure. However, the effects of 
changing r are covered, to some extent, by those of e. The impact parameter b is not essential, 
as seen in the following expression for the total emission power W of SPR, whose b dependence 
is collected into a simple scaling law [19] 




(6) 



Pem{CO,k z )\ b = e- 2 ^ b - h ^P em {(0,k z )\ bo , (7) 

where P em (a>,k z )\i, is the (0- and ^--resolved emission power for an impact parameter b and 
bo is a reference impact parameter chosen arbitrarily. Therefore, uSPR and cSPR change in a 
straightforward way as b varies, with the underlying physics unaltered. In the following sub- 
sections, therefore, five parameters, v,e,N,N[, and 0, are varied in this order to see how each 
affects the spectrum. 

4.1. Velocity 

The velocity of the electron beam is a key parameter in the uSPR. Indeed, as seen in the previous 
section, the SPR at v = 0.5c is understood using the theory of cSPR, while at v = 0.99999c the 
uSPR also plays an important role. We shall examine how the conventional picture fails with 
varying electron velocity. An obvious but nonessential v-dependence is an increase of the SPR 
intensity due to the v dependence of the decay-constant |T|; if impact parameter b is fixed, the 
overall SPR spectrum behaves as exp(— 2|r|£>). To eliminate this trivial v-dependence, we have 
set 

v 1 

b = 0MBva, with B = -, y = , (8) 

considering |T| oc l/(/3y) for k z = 0. 

The reflected intensity maps for the monolayer of = 21 are shown in Figs. 9 (a) and (b) 
for v = 0.7c (y = 1 .4) and in (c) and (d) for 0.99c (y = 7.09), along with the PhB structure. 
Panels (a) and (c) show only the SPR intensity, while they are superposed by the PhB structure 
in panels (b) and (d). 

At v = 0.7c, there is a marked bright line along the shifted v line of h = 1 . Along the line, the 
intensity contrast of the SPR is quite strong at low frequencies. In particular a point-like reso- 
nance is seen at coa/2nc ~ 0.745. As panel (b) shows, this resonance arises just at a crossing 
between the dispersion curve of an even-parity PhB and the shifted v line of h = 1 . Therefore, 
this is a type (A) signal according to the classification of the last section. 

We can see the flat streaks of strong intensity just at this frequency,. We note that the signal 
becomes stronger as k x approaches the backward light line co = —ck x . This feature is common 
to all the horizontal streaks appearing at the frequencies of the pseudo gaps. These are signals 
of type (E) of the uSPR. We should note that the PhB mode, which crosses the shifted v lines, 
has a negative group velocity, and the excited PhB mode propagates in the — x direction. A 
backward-oriented diffraction taking place at the left edge of the PhC explains the tendency 
towards the line (0 = —ck x . Analogous flat lines exist, for instance, at (Oa/lnc ~ 1.09. 

In addition, we see clearly a high-intensity spectrum appearing almost parallel to the shifted 
v lines. The curves are in fact coincident with the dispersion curves of quasi-guided PhB of 
the even parity. Therefore, they are type (B) signals. Note that the odd-parity PhB dispersion 
curves are also visible, with reduced strength as compared to the even-parity PhBs. Altogether, 
at v = 0.7c, cSPR coexists with uSPR and odd-parity PhBs are seen in the uSPR spectrum, 
with weaker intensity than even-parity PhBs, however. Combining this result with what we 
have seen in Sec. Ill for v = 0.5c and v = 0.99999c, we may conclude that, as v increases from 
v = 0.5c, the uSPR becomes visible and the even-parity selection rule of uSPR is less stringent 
at non-ultra-relativistic velocities. 

The result for v = 0.99c indeed confirms this conclusion. At v = 0.99c, several bright curves 
arise in Figs. 9 (c) and (d) with little intensity contrast along the PhB dispersion curves. This is 
the type (B) signal. We can observe odd-parity excitation of weak intensity. Therefore, although 
the even-parity selection rule is indeed dominant, it is somewhat relaxed for v = 0.99c. On the 
shifted v line, there are signals of cSPR, as theory predicted for type (A) features in Sec. III. 




Fig. 9. Dependence of the reflected SPR spectra on electron velocity. The result for a finite 
monolayer (N = 21) of contact cylinders is shown for e = 2.05. Panels (a) and (b) are the 
results for v = 0.7c, and panels (c) and (d) are those for v = 0.99c. Panels (b) and (d) are 
reproductions of (a) and (c), respectively, overlaid with the PhB structure of N = °°. The 
PhB modes with even (odd) parity are indicated by red (green) circles. See text for the 
impact parameter b used in the calculation. 



The breakdown of the even-parity selection rule at these velocities is explained as follows. 
With a decrease of v, |T| increases to make the incident evanescent light decay more quickly 
when passing the monolayer. This increases the asymmetry of the evanescent wave with respect 
to the mirror plane and makes the even-parity selection rule less effective. The degree of the 
symmetry of the input wave may be given by the factor exp( — |T|2r), called here the symmetry 
factor, which measures the decay of the evanescent wave while traversing the PhC in the +y 
direction. If this factor is unity, the evanescent light seen by the PhC is mirror-symmetric. At 
v = 0.99c, the symmetry factor is 0.408 at coa/2nc = 1 and too small to guarantee strictly the 
even-parity selection rule. Therefore, odd-parity PhBs are allowed somewhat as uSPR signals. 

The results for the other v are briefly summarized without giving numerical results. At v = 
0.9c, when the symmetry factor is 0.047, cSPR and uSPR coexist and odd-parity PhBs are 
seen in the latter. At v = 0.999c, the asymmetry factor increases to 0.755. The intensity map 
gradually tends to the case of v = 0.99999c with the symmetry factor 0.972; signals along the 
odd-parity PhBs disappear, leaving behind only the even-parity signals as type (B) signals. The 
horizontal bright streaks appear solely in the regions of the pseudo-gaps of even-parity bands. 

Finally, we should comment on the case of non-zero v z . The critical velocity of the electron, 
above which the uSPR begins to emerge does not change so much by non-zero v,. An important 
point is that at ultra-relativistic velocities the evanescent wave can be effectively regarded as a 
plane wave. This is not controlled by v z , but is controlled by v, the magnitude of the velocity 
vector. However, other features of the uSPR changes as discussed in the previous section. 

In conclusion, in the frequency region coa/2nc ~ 1, the uSPR is conspicuous when v exceeds 



0.7c, and the even-parity selection rule holds progressively better as v approaches c from 0.9c. 

4.2. Dielectric constant 

For a PhC with r and N kept fixed at r = 0.5a and N = 21, let us examine how the SPR 
spectrum varies as the dielectric constant e of the cylinders changes in the monolayer. We 
select three values of e, e = 4.41, 1 — ((O p /(o) 2 and — °°. The first case corresponds to the 
dielectric constant of fused quartz with £ nearly twice as large as that used above, the second 
is the dielectric constant of a Drude metal with a> p the plasma frequency, and the third is the 
dielectric constant of a perfect conductor. To avoid the poor convergence of the cylindrical- 
wave expansion for the metallic cylinders in contact, we created a narrow opening between the 
cylinders by setting r = 0.45a in the Drude case. We assumed co p a/2nc = 1, i.e., the plasma 
wavelength equals the lattice constant. Calculation is made for the monolayer system using 
v = 0.99999c and b = 3.33a, as before. 

The reflected SPR intensity maps are shown in Fig. 10, together with the corresponding 
PhB structure obtained for the N = °° system. Panels (a) and (b) depict the result of dielectric 
cylinders, panels (c) and (d) treat the Drude cylinders, and panel (e) presents the spectrum of 
the cylinders of a perfect conductor. Panels (b) and (d) also involve the band structures of the 
monolayer. Considering v is ultra-relativistic, we only plotted the even-parity PhB structure. 
Note that for the perfect conductor case of panel (e), the ODOS does not present any peaks 
except for Wood's anomaly and the PhB structure is completely absent. 

Clearly, the calculated uSPR intensity shown in panels (a) and (b) is well correlated with the 
PhB structure. Namely, the bright curves of strong SPR intensity are type (B) signals having a 
positive slope and tracing very well the even-parity PhB dispersion curves. In addition, we can 
recognize type (E) signals of the bright flat lines terminated at the backward light line, which 
are seen just at the frequencies of the pseudo-gaps. These features agree with what we have 
seen in Sec. III. Most importantly, the PhB structure with larger dielectric constant is indeed 
probed by the uSPR. 

As for the Drude case shown in panels (c) and (d), the PhB structure is composed of many 
flat bands. These PhB modes have their origin in the tight-binding coupling among cylinders 
of the surface plasmon polaritons (SPP), localized on each cylinder surface [29, 30, 31]. The 
calculated intensity map demonstrates that these SPP bands are coupled only weakly to the 
incoming evanescent wave. In contrast, the PhB around (Oa/lnc = 0.5, which has a modest 
group velocity, is strongly coupled to the evanescent wave, yielding a very strong SPR signal. 
We thus conclude that uSPR carries information of PhBs of SPP origin. 

Finally, in Fig. 10(e) the strong intensity of the SPR arises solely on the shifted v lines. This 
reflects the absence of the PhB structure in the array of perfect-conductor cylinders. Thus, we 
can conclude that the uSPR is peculiar to dielectric and metallic PhCs with finite dielectric 
function and is completely absent in the systems without PhBs. 

4.3. Number of cylinders 

In the numerical results shown so far, the number of the cylinders is fixed at either N = 21 or 
N = oo. At small N, typically less than 8, the PhB structure is not clearly visible in the intensity 
map of the SPR on the (k x , Co) plane. On the other hand, at large N the PhB structure is clearly 
visible as demonstrated in Figs. 5 and 6. In this region, however, change of the SPR intensity 
map with increasing N is less remarkable. Nevertheless, if we have a close look at the spectral 
line shapes of the uSPR, they indeed change with N. To investigate this feature, we consider 
the SPR spectra at a fixed solid angle (0,(j>) = (60°, 180°) as a function of frequency. Figure 
1 1 shows the spectra for various N. The SPR signals are strongly enhanced around a>a/2nc = 
1.195, which corresponds to an intersection point between the line of k x —h= (o/c)cos0 





Fig. 10. Reflected SPRs from the monolayers of contact cylinders of various dielectric 
constants. Panel (a) shows the result of dielectric cylinders of e = 4.41 and r = 0.5a, panel 
(c) is the result of metallic cylinders of Drude dielectric constant e = 1 — a>~/(0 with 
(OpCi/lnc = 1, and panel (e) shows the result of cylinders of perfect metal, i.e., r = 0.5a 
and e = — oo. Panels (a) and (c) are reproduced in panels (b) and (d), respectively, with the 
corresponding PhB superposed. The same parameters as in Fig. 4 were used for the electron 
beam. 



30 




(Oa/27tc 

Fig. 11. The reflected SPR intensity spectra at (d,(j>) = (60°, 180°) from finite-size PhCs 
of various N. The same parameters as in Fig. 4 were used except for N. 

(see Eq. (3)) and the PhB structure co = (O n (k x ,0). The intersection point is off the shifted 
v lines and thus is indeed an uSPR signal. We can clearly observe that as N increases, the 
intensity at the peak grows but seems to be saturated to a certain value. This implies that the 
radiation intensity of the uSPR per unit length of the electron trajectory decreases at large N 
and eventually vanishes at N = °°. This property is reasonable because the uSPR is completely 
forbidden in the system of perfect translational invariance with N = °°. On the other hand, 
we also found that the intensity of the cSPR on the shifted v lines increases almost linearly 
with N, as expected from the conventional theory of the SPR. Therefore, at very large N the 
cSPR signals will dominate over the uSPR ones. However, even at N = 200 we found that 
the SPR intensity map does not differ so much from Fig. 5, in which the uSPR signals are 
rather stronger than the cSPR ones. Besides, in Fig. 1 1 we can clearly observe that the spectral 
width of the peak decreases with increasing N. This property reflects better confinement of the 
radiation energy for larger N. This width should converge to a certain value atN = °°, which is 
inversely proportional to the life-time of the relevant photonic band mode. This is nothing but 
the homogeneous broadening of the spectral line width of the uSPR. 

4.4. Number of stacking layers 

So far we have confined ourselves to the monolayer PhC. Now let us stack the identical mono- 
layers periodically in the y direction. As we increase the number N/ of the stacking layers, the 
ODOS reveals a progressively finer structure as a function of frequency. Each peak of ODOS 
corresponds to a quasi-guided PhB mode confined in the stacked layers. The typical peak-to- 
peak distance in frequency is inversely proportional to N/. Moreover, each peak is generally 
getting sharper. 

If N/ is large enough, the scattering of the evanescent wave in the ultra-relativistic regime is 
identical to the transmission and reflection of a TE-polarized plane wave that enters the PhC 
with its left edge as an entrance surface. The slab PhC in question has a finite thickness Na in 
the x direction and has a large extension in the y direction with the entrance surface parallel to 
the yz plane. The wave vector component k\\ parallel to the entrance surface is conserved, and 
the incident plane wave excites the bulk PhB modes having the same k\\ . There is no momentum 
conservation in the x component; in principle the light excites any PhB modes of arbitrary k x . 
The scattered wave is decomposed into diffraction channels [27] . 




Fig. 12. Reflected SPR intensity map of various PhCs of cylinders. Panel (a) shows the 
result of the double-layer system of contact low-index cylinders (N = 21 and N/ = 2). 
Panel (b) is produced from panel (a) by overlaying the even-parity PhBs. The mirror plane 
of the parity lies in between the double-layers. Panel (c) shows the result for a multi-layer 
PhC of square lattice of contact low-index cylinders (N = 8 and N; = 20). Panel (d) is a 
reproduction of panel (c), overlaid with the PhB structure along the T — X direction. Only 
the even-parity modes with respect to the mirror plane relevant to T — X are shown. The 
same parameters as in Fig. 4 were used for the electron beam. 

Let us examine the stacked monolayers with a square lattice of the cylinders. For uSPR, 
(±F,k z ) plays the role of k\\ in the above analogy, and thus we may set k\\ — in the ultra- 
relativistic case, provided k z = 0. Accordingly, the incident plane wave has the wave vectors 
(k x , 0, 0) and propagates in the r — X direction of the square lattice. The bulk PhB modes along 
r — X are thus excited. The conclusion is thus that uSPR will carry information along r — X of 
the bulk PhB modes if both N and N/ are sufficiently large. According to the above arguments, 
the SPR intensity is expected to be enhanced in the forward and backward directions, which 
correspond to the specular transmission and reflection. In addition, if coa/2nc > 1, the SPR 
intensity is expected to be enhanced also on the curves 




which correspond to the diffraction channels associated with the reciprocal lattice points h = 
2nn/a (n : integer) along the k y axis. 

Figure 12 illustrates the reflected SPR intensity maps and the relevant PhB structure of low- 
index cylinders in contact. We used e = 2.05, v = 0.99999c and b = 3.33a. In Fig. 12(a) the 
spectrum from the double-layer (Ni — 2) structure with N = 21 is shown. The intensity map 
overlaid with the PhB structure of the double layers (but for = °°) is shown in panel (b). As 
before, we plotted only the even-parity PhB structure. In the double layer, the mirror plane lies 



midway between the layers. We see the number of bands is almost twice that of the monolayer 
band structure shown in Fig. 3. This is reasonable, since the degenerate band-structures of each 
of the monolayers are split in the double layer. Obviously there is a very good correlation of the 
strong signals of uSPR with the band structure of the even parity. 

Figure 12(c) shows the reflected intensity map of the finite multilayers PhC of N/ = 20 and 
N = 8. We consider this to be a test system simulating the slab-type PhC of square lattices. 
We observe at once a signal of high intensity along a hyperbolic curve whose bottom is found 
at (k x a/2n, coa/2nc) = (0, 1). Obviously, this curve corresponds to Eq. (9) with /; = 1. Strong 
SPR signals other than the hyperbolic curve are found at (Oa/lnc = 0.73, 0.93, and 1.46. To 
identify these signals, Fig. 12(c) was overlaid with the even-parity PhB structure along the 
r — X direction of the square lattice. The result is shown in Fig. 12(d). As can be clearly seen, 
the strong signals correspond to the anti-crossing points of the even-parity PhB structure. The 
bright curve connected to the strong signal around coa/2nc = 0.73 is shown to be along the 
PhB dispersion curve. Thus, we can conclude that the intensity map of the uSPR correlates 
well with the corresponding PhB structure even in the case of stacked monolayers. 

4.5. Azimuthal angle 

So far, we have considered the case of k z = (<p = 0° and 180°), that is, we have examined the 
radiation emitted within the xy plane. We here investigate the dependence. For this purpose, 
we write the differential cross section of SPR in polar coordinates [19] 

dw _?V?^ (|/ m (0O|2 + |/> , )|2)j (10) 
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Obviously, the cross section must have inversion symmetry under the operation <j> — > — 0, re- 
flecting the inversion symmetry of the z coordinate with respect to the electron trajectory located 

atz = z - 

The radiation emission of non-vanishing lc z is generally small compared with that of k z = 0, 
and k z is a conserved quantity in the scattering by the PhC. Therefore, the k z dependence of the 
observed SPR will be controlled dominantly by that of the decaying exponential exp(— \T\b) of 
the initial light. This exponential decreases with increasing \k z \, so that the radiation is domi- 
nated by the SPR of k z = 0. 

Let us consider the radiation emission toward solid angle {9,(p). In the ultra-relativistic 
regime it follows that 



co\ 2 fC0\ 2 ,co . 

— — ( — J — (— sin0sin0)- ~ iqsinO sin0 . (13) 
c ) \ v / c 

Thus, the SPR cross section at a given co and 0°, 180°) is dominated in the forward (8 = 0°) 
and backward (9 = 180°) directions. Similarly, at a given 9, the SPR cross section is dominated 
around the plane perpendicular ((j) ~ 0° and 180°) to the cylindrical axis. 

Figure 13(a) projects the reflected SPR spectrum Eq. (10) for <p = 165° onto the (k x , co) plane, 
and Fig. 13 (b) projects the reflected SPR spectrum for 9 = 90° is projected onto the (k z , co) 
plane. Figure 13(a) verifies that the strong intensity is limited around the forward and backward 
light lines, as asserted above. The intensity contrast along the backward light line is related to 
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Fig. 13. Reflected SPR intensity off the xy plane. Spectra from the monolayer of contact 
low-index cylinders are shown, (a) Intensity in the (k x ,co) plane at <j> = 15°. (b) Intensity 
in the (k~,co) plane at 6 = 90°. For the electron beam, the same parameters are used as in 
Fig. 4. 



the PhB structure with finite k z . In Fig. 13(b), the strong SPR intensity is limited around k z = 0. 
At k- = 0, three marked peaks can be found at CO — 0.75, 1, and 1.5. They correspond to the 
crossing points between the bright curves of Fig. 5 and the line of k x = (i.e., 6 = 90°). Thus, 
we can conclude that SPR is highly directive within the plane normal to the cylindrical axis. 

5. Summary and discussions 

To summarize, we have presented a theory of uSPR that arises when an ultra-relativistic electron 
beam is used to obtain the SPR from a finite PhC composed of cylinders. The ultra-relativistic 
electron accompanies an evanescent wave whose spatial decay is almost negligible at k z = 0, 
so that the evanescent wave can be regarded as a plane wave propagating in the direction of the 
trajectory. This yields a peculiar radiation emission from the PhCs, which cannot be explained 
by the conventional theory of the SPR in which the finiteness of PhC is treated as infinite. The 
spectrum of the uSPR can be used as a probe of the PhB structure of the quasi-guided modes 
having the even-parity symmetry with respect to the relevant mirror plane. 

We have also presented properties of the uSPR in detail by changing several system parame- 
ters. We found that the uSPR coexists with the cSPR at moderate velocities typically in between 
0.7c and 0.99c. We also found that the uSPR is completely absent in the perfect-conductor 
cylinders because of the absence of PhBs. Otherwise, the spectra of the uSPR correlate with the 
corresponding PhB structure very well. We also found that the cross section of the SPR at an 
ultra-relativistic velocity is highly directive within the plane normal to the cylindrical axis. 

It should be emphasized again that the uSPR is an unexpected phenomenon, in which the 
conventional theory assuming infinite periodicity of the PhCs fails to reproduce its features. On 
the other hand, the present theory explains very well both the c and u SPRs in a unified manner. 
There are three key items in the uSPR: presence of PhBs, broken translational invariance of 
PhC, and ultra-relativistic velocity of the electron beam. Lack of either one of the three items 
prevents understanding of the uSPR correctly. 

In actual PhCs various types of disorder or randomness are inevitable, yielding the inho- 
mogeneous broadening of spectral line width of the SPR. For instance, a rigid vibration of 
constituent cylinders of the PhC gives rise to the Brillouin scattering. As a result, the broad- 
ening of the line width is given by the frequency of the vibration. The relative percentages of 
various disorder factors depend crucially on the frequency range concerned. Thus, when we 
extract the intrinsic SPR signals from PhC, we should carefully take account of disorder. 



From the point of view of the radiation emission from high-energy electron interacting with 
periodic structure, we should comment on the peculiarity of the uSPR in comparison to the 
channeling radiation, or in other words, Kumakhov radiation [32]. The latter radiation occurs 
inside a crystal when a certain condition is satisfied for the incident angle of the electron with 
respect to a major crystal direction. The radiation depends strongly on the meandering tra- 
jectory of the electron trapped around a crystal plane or a crystal axis and has a monotonic 
frequency. On the other hand, the uSPR does not require the meandering of the electron tra- 
jectory. Actually, in our theoretical approach, the trajectory of the electron is assumed to be 
straight. In addition, the radiation spectrum of the uSPR is not monotonic for a fixed trajectory 
and the typical frequency range is inversely proportional to the lattice constant of the PhC under 
consideration. Therefore, the uSPR is not categorized into the channeling radiation. 

The uSPR is, in some sense, similar to the transition radiation [33] because the broken trans- 
lational invariance along the electron trajectory is crucial in both the radiations. However, there 
is a marked difference in the directivity between the transition radiation and the uSPR. Suppose 
that an ultra-relativistic electron passes from vacuum to a dielectric medium. It is better to focus 
on the radiation into the vacuum side, because the induced radiation in the medium is a mixture 
of the transition radiation and the Cherenkov radiation. It was shown that this radiation into the 
vacuum side is backward-oriented. On the other hand, as we showed in the paper, such a high 
directivity into the backward direction is only possible if the relevant frequency lies in a pseudo 
gap of the PhB structure. Among other electron-induced radiations, the uSPR may have the 
closest resemblance to the diffraction radiation regarding the broken translational invariance 
and the trajectory which does not pass through any air/dielectric interfaces. To further clarify 
the resemblance, a detailed investigation of the diffraction radiation in PhC is in order. 
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